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Magnetic and structural properties of ferromagnetic Fe5PB2 and Fe5SiB2 and effects of 

Co and Mn substitutions 


Michael A. McGuire* and David S. Parker 
Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831 USA 

Crystallographic and magnetic properties of Fe5PB2, Fe4CoPB2, Fe4MnPB2, Fe5SiB2, Fe4CoSiB2, 
and Fe4MnSiB2 are reported. All adopt the tetragonal CrsBa structure-type and are ferromagnetic 
at room temperature with easy axis of magnetization along the c-axis. The spin reorientation 
in Fe6SiB2 is observed as an anomaly in the magnetization near 170 K, and is suppressed by 
substitution of Co or Mn for Fe. The silicides are found to generally have larger magnetic moments 
than the phosphides, but the data suggests smaller magnetic anisotropy in the silicides. Cobalt 
substitution reduces the Curie temperatures by more than 100 K and ordered magnetic moments by 
16 - 20 %, while manganese substitution has a much smaller effect. This suggests Mn moments align 
ferromagnetically with the Fe and that Co does not have an ordered moment in these structures. 
Anisotropic thermal expansion is observed in Fe5PB2 and Fe5SiB2, with negative thermal expansion 
seen along the c-axis of Fe5SiB2. First principles calculations of the magnetic properties of Fe5SiB2 
and Fe4MnSiB2 are reported. The results, including the magnetic moment and anisotropy, and are 
in good agreement with experiment. 


I. INTRODUCTION 

Uniaxial ferromagnets composed of earth-abundant el¬ 
ements are of interest to researchers exploring new ma¬ 
terials for permanent magnet applications [1]. Fe 5 PB 2 
and Fe 5 SiB 2 are examples of such materials that have 
not been heavily investigated in this respect to date. 
The crystal structures and compositions of these com¬ 
pounds were first investigated by Aronsson and cowork¬ 
ers [2, 3] and Rundqvist [4]. They were shown to be 
isostructural, adopting the tetragonal CrsBa structure 
type. Ferromagnetism was first reported in the phosphide 
based on magnetization measurements [5] and Mossbauer 
spectroscopy [6], which identified the uniaxial nature of 
the magnetic structure with the easy axis of magnetiza¬ 
tion along the crystallographic c-axis. Soon thereafter, 
Fe 5 SiB 2 was identified as a uniaxial ferromagnet at room 
temperature as well [7], and a spin reorientation was 
seen to occur near 140 K, resulting in planar magnetic 
anisotropy at lower temperature [8]. The Curie temper¬ 
atures (Tc) of both phases were found to be well above 
room temperature, near 630 K for Fe 5 PB 2 [5] and 784 
K for Fe5SiB2[7]. A very recent study of Fe 5 PB 2 reports 
Tc = 655 K, magnetic anisotropy comparable to hard 
ferrites, and a magnetic moment per Fe of 1.72 hb [9]. 
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The isostructural Mn compounds are also ferromagnetic, 
but with Curie temperatures closer to room temperature, 
making them perhaps more interesting as magnetocaloric 
rather than magnet materials [10, 11]. No reports of 
magnetization measurements for Fe 5 SiB 2 were located in 
the literature, and no detailed crystallographic study has 
been published to date. In addition, understanding of 
effects of chemical substitutions in these uniaxial, high- 
temperature ferromagnets is presently lacking. 


In this work, we report structural and magnetic prop¬ 
erties determined from measurements on polycrystalline 
samples of Fe 5 PB 2 , Fe 5 SiB 2 , Fe4CoPB2, Fe4CoSiB2, 
Fe 4 MnPB 2 , and Fe 4 MnSiB 2 . Full crystal structure re¬ 
finements for all compounds are reported at room tem¬ 
perature, and for both Fe 5 PB 2 and Fe 5 SiB 2 at 20 K. Sig¬ 
nificant anisotropy is observed in the thermal expansion, 
which is negative along the c-direction in Fe 5 SiB 2 . Mag¬ 
netization measurements are used to determine the sat¬ 
uration magnetizations and identify the Curie and spin- 
reorientation temperatures. Results of density functional 
theory calculation of the ordered magnetic moment and 
anisotropy of Fe 5 SiB 2 and Fe 4 MnSiB 2 are reported, and 
compared to the measurements. The silicides are found 
to have a larger magnetic moment than the phosphides, 
but smaller magnetic anisotropy at 300 K is indicated in 
the silicides. Both chemical substitutions suppress the 
spin reorientation in Fe 5 SiB 2 , indicating Fe 4 CoSiB 2 and 
Fe 4 MnSiB 2 are uniaxial ferromagnets down to at least 
30 K. In both Fe 5 PB 2 and Fe 5 SiB 2 , substitution of Mn 
expands the unit cell, has little effect on the Curie tem¬ 
perature, and produces a slight reduction in magnetic 
moment. Substitution of Co contracts the unit cell, and 
strongly suppresses both the Curie temperature and the 
magnetic moment. These findings suggest that the or¬ 
dered magnetic moment on Co in these structures is ei¬ 
ther very small or aligned non-collinearly with the net 
magnetization. 
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II. EXPERIMENTAL PROCEDURES 

Polycrystalline samples of nominal compositions listed 
in Table I were synthesized by reacting the elements in 
evacuated and flame-sealed silica tubes. The tubes were 
heated to 1000 °C at 100 °C per hour, and held at this 
temperature for 12 hours. The resulting mixtures were 
then ground into powders and pressed into 1/2 inch di¬ 
ameter pellets at room temperature. The pellets were 
sealed in evacuated silica tubes and heated at 1000 °C 
for one week. The pellets were then ground into pow¬ 
der, re-pelletized, and heated in evacuated silica tubes 
at 1000°C for 5 days. Samples produced in this way 
were found to be ^90% pure. For the silicon-containing 
samples (with no phosphorus), the initial mixing of the 
elements can be done by arc-melting, but subsequent 
grinding, pelletizing, and annealing are needed to im¬ 
prove the phase purity. This later route was used to 
make the Fe 4 CoSiB 2 used in this study, because it gave 
a final sample with less ferromagnetic impurity. 

X-ray diffraction was performed on a PANalytical X- 
Pert Pro diffractometer using monochromatic Cu Kai 
radiation in Bragg-Brentano reflection geometry. An 
Oxford PheniX cryostat was used for diffraction mea¬ 
surements below room temperature. Quantitative anal¬ 
ysis and crystal structure refinement was done with Ri- 
etveld refinement using the program FullProf [12]. To 
determine the easy axis of magnetization by diffrac¬ 
tion, finely ground samples were sprinkled onto a vac¬ 
uum grease coated glass slide with a strong perma¬ 
nent magnet underneath. Compositions of the primary 
phases were examined using energy dispersive analysis 
with a Bruker Quantax 70 x-ray spectrometer (EDS) 
in a Hitachi TM3000 tabletop scanning electron micro¬ 
scope. Magnetization measurements on solid polycrys- 
talline samples cut from the sintered pellets were per¬ 
formed with a Quantum Design Magnetic Property Mea¬ 
surement System. 

III. RESULTS AND DISCUSSION 
A. X-ray diffraction 

Analysis of powder x-ray diffraction patterns (Fig. 1) 
shows that all of the materials studied adopt the tetrago¬ 
nal CrsBa structure type at room temperature. The sam¬ 
ples used in the current study were ^90% pure, and the 
impurities identified in the diffraction patterns along with 
their weight percent concentrations are shown on the Fig¬ 
ure. In the Fe 5 PB 2 -based materials, the only secondary 
phase detected was Fe 2 P, which is paramagnetic at room 
temperature, and undergoes a first-order ferromagnetic 
phase transition near 216 K [13]. About 3 % of Fe 5 _ 5 Si 2 B 
[2, 3] was observed in the Fe 5 SiB 2 and Fe 4 MnSiB 2 sam¬ 
ples. The magnetic properties of this phase have not 
been reported, but in the course of this study we have 
synthesized and performed preliminary characterization 


of a sample of composition Fe 4 . 7 Si 2 B and find that it is 
likely paramagnetic at room temperature and antiferro¬ 
magnetic at lower temperature. The Fe 4 CoSiB 2 sample 
(Fig. le) contains a small amount of an unidentified im¬ 
purity phase, which may be reflected in the high temper¬ 
ature magnetization data presented below. 

Crystallographic data determined from refinement of 
the diffraction data are collected in Table I. The struc¬ 
ture determined for Fe 5 PB 2 and listed in Table I is 
in good agreement with that originally reported by 
Rundqvist and coworkers [4, 6 ] and confirmed recently 
with single crystal x-ray diffraction [9]. The unit cell pa¬ 
rameters of Fe 5 SiB 2 in Table I are also consistent with 
literature reports [3], although a full structure refinement 
has not been previously reported. The crystal structure 
is shown in Fig. 2b. Sites labeled Ml and M2 are oc¬ 
cupied by the transition metal atoms. At z = 0 and ^ 
are planes of composition M^B 2 . Slabs of composition 
MI 4 P fill the space between these planes. The Ml posi¬ 
tions has site symmetry m and an irregular coordination 
environment. It is coordinated by two P atoms, three 
B atoms, two M2 sites, and two Ml sites. M2 with 
site symmetry 4/m is at the center of a square of four 
B atoms in the 06 -plane and a slightly distorted cube of 
eight Ml atoms. 

The ratio of transition metals in the main phase of 
the Co and Mn substituted samples was determined by 
EDS and is shown in Table I. This was determined by 
averaging multiple measurements on single grains of the 
primary phase. The results are consistent with the nom¬ 
inal ratio of 1:4. Allowing mixing of P or Si and B on 
their respective 4o and 86 sites did not result in any im¬ 
provement of the fits to the diffraction data for any of 
the samples studied here, but it did indicate that some 
B may replace the heavier main-group elements, as re¬ 
ported previously in Fe 5 PB 2 [ 6 , 9]. For the final struc¬ 
tural rehnements occupancies were fixed at the nominal 
composition. 

Substitution of Co is seen to reduce the lattice param¬ 
eters and unit cell volume, while substitution of Mn has 
the opposite effect. The c/a ratio is not substantially af¬ 
fected by the substitutions, but it is slightly different be¬ 
tween the P- and Si-containing materials. These changes 
are as expected based on the published lattice constants 
of C 05 PB 2 (5.42 A, 10.20 A) [4], Mn 5 PB 2 (5.54 A, 10.49 
A) [ 4 ], and Mn 5 SiB 2 (5.61 A, 10.44 A) [2]. Note that 
Co 5 SiB 2 is not known to form. 

All of the materials are strongly ferromagnetic at room 
temperature {vide infra). Mossbauer spectroscopy has 
shown that Fe 5 PB 2 and Fe 5 SiB 2 are uniaxial ferromag- 
nets at room temperature, with the easy axis of mag¬ 
netization along the tetragonal c-axis [ 6 - 8 ]. To confirm 
this in the samples studied here, and to determine if the 
alloying of Co or Mn may change the easy axis, diffrac¬ 
tion patterns were collected from powders aligned in a 
magnetic held so that the easy axis of magnetization was 
normal to the sample plane. The results are shown in 
Fig. 2a. Three rehections are observed in the angular 
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FIG. 1. Room temperature powder x-ray diffraction patterns and fits from Rietveld refinements. Difference curves are shown 
at the bottom of each panel. The lower set of tics in panels a-d and f locate Bragg reflections from the impurity phases noted on 
each plot. The impurity in Fe 4 CoSiB 2 has not been identified and its concentration is estimated from relative peak intensities. 


TABLE I. Crystallographic and magnetic properties. All materials adopt the CrsBs structure type with space group 74/mcm 
and Ml at 16( (x, x-|-l/2, z), M2 at 4c (0, 0, 0), P/Si at 4a (0, 0, 1/4), B at 8h (x, 1/2-l-x, 0). Note J = fioM, and divide J in 
Tesla by dvr x to convert to M in emu/cm®, or by 4^ x 10“^ to convert to M in A/m. 


nominal composition 

Fe5PB2 Fe4CoPB2 Fe4MnPB2 

FesSiBa Fe4CoSiB2 Fe4MnSiB2 

EDS (Co,Mn)/Fe ratio 

- 

0.28 

0.24 

- 

0.25 

0.26 

a (A) 

5.4867(1) 

5.4522(1) 

5.4947(1) 

5.5507(1) 

5.5332(2) 

5.5689(1) 

c (A) 

10.3532(2) 

10.3585(3) 

10.3754(2) 

10.3359(2) 

10.2692(4) 

10.3582(2) 

V (A®) 

311.67(1) 

307.92(1) 

313.25(1) 

318.45(1) 

314.40(2) 

321.24(1) 

density (g / cm®) 

7.075 

7.227 

7.020 

6.864 

7.017 

6.785 

c/a 

1.887 

1.900 

1.888 

1.862 

1.856 

1.860 

xmi 

0.1701(4) 

0.1698(5) 

0.1699(4) 

0.1678(4) 

0.1694(5) 

0.1685(4) 

zmi 

0.1400(3) 

0.1419(4) 

0.1410(3) 

0.1384(3) 

0.1387(4) 

0.1381(3) 

Xb 

0.384(4) 

0.385(5) 

0.388(4) 

0.386(4) 

0.372(5) 

0.382(4) 


1.2 

1.3 

1.5 

1.3 

1.2 

1.4 

Js 5 K (T) 

1.20 

1.03 

1.18 

1.31 

1.06 

1.22 

Ms 5 K (Ats/F.U.) 

8.00 

6.71 

7.51 

9.06 

7.24 

8.52 

Js 300 K (T) 

1.09 

0.87 

1.00 

1.21 

0.961 

1.17 

Ms 300 K (/tb/F.U.) 

7.29 

5.68 

6.72 

8.39 

6.55 

8.14 

Tc (K) 

640 

515 

650 

> 780 

675 

770 
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FIG. 2. (a) Powder x-ray diffraction results showing the re¬ 
sults of alignment of the powders with a magnetic field normal 
to the sample surface, (b) The crystal structure adopted by 
the materials studied here. 


range shown. They are labeled by their indices on the 
top panel. For each material, the intensity of the 004 re¬ 
flection is strongly increased relative to the 022 and 121 
reflections when the powders are magnetically aligned. 
This indicates that the easy axis of magnetization is 
parallel to the c* reciprocal lattice vector. For tetrag¬ 
onal symmetry, this is the c-axis of the crystal structure. 
Some intensity is still observed on the 022 and 121 re¬ 
flections of the aligned samples, which is expected if the 
alignment is incomplete or if all of the powder particles 
are not single crystals, but instead are composed of ag¬ 
glomerates of misaligned crystallites. 

Powder x-ray diffraction patterns from Fe5PB2 and 
Fe5SiB2 were also collected at 20 K. The results from 
both materials were well described by the room tem¬ 
perature structures, with lattice parameters 5.4766(1) A 
and 10.3488(2) A for the phosphide and 5.5412(1) A and 
10.3390(2) for the silicide. Other crystallographic param¬ 
eters determined at 20 K are shown in Table II. A crys¬ 
tallographic distortion may be expected to occur at the 



r(K) 


FIG. 3. Temperature dependence of (a) the positions of 
the 130 and 006 Bragg reflections from Fe 5 PB 2 and Fe 5 SiB 2 
and (b) the magnetization measured approximately along 
the room temperature easy axis for Fe 5 SiB 2 , Fe 4 GoSiB 2 , 
Fe4MnSiB2, and Fe5PB2. 


spin reorientation temperature in Fe5SiB2; however, no 
distortion was resolved in the data, indicating the magne¬ 
toelastic coupling in the a6-plane is not strong. Compar¬ 
ison of lattice constants at room temperature (Table I) 
and 20 K (Table II) reveal strongly anisotropic behavior. 
In Fe5PB2, the contraction of a upon cooling from room 
temperature to 20 K is more than four times larger than 
the relative contraction along c. In Fe5SiB2, the basal 
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TABLE II. Crystallographic properties at 20 K in the CrsBs 
structure type with space group /4/mcm and Fel at 161 (x, 
x+1/2, z), Fe2 at 4c (0, 0, 0), P/Si at 4a (0, 0, 1/4), B at 8 h 
(x, 1/2+x, 0). 


Fe5PB2 Fe5SiB2 


a 

5.4766(1) 

5.4812(1) 

c 

10.3488(2) 

10.3390(2) 

V 

310.39(1) 

317.46(1) 

XFel 

0.1698(3) 

0.1676(3) 

ZFel 

0.1407(2) 

0.1379(2) 

Xb 

0.388(2) 

0.388(3) 


1.8 

1.5 


plane also changes more strongly with temperature than 
the c-axis. In fact, the c-axis shows negative thermal ex¬ 
pansion in this compound. The thermal expansion was 
examined further by monitoring the 130 and 006 Bragg 
reflections upon cooling from 300 to 20 K. The results are 
shown in Figure 3a. The position of the 130 reflection de¬ 
pends only on the basal plane lattice constant a and the 
006 position depends only on the c-axis length. Normal 
behavior is seen in Fe5PB2. The negative thermal expan¬ 
sion along c is apparent in the temperature dependence of 
the position of the 006 reflection in Fe5SiB2. This com¬ 
pound is known to undergo a spin reorientation below 
room temperature, with the easy axis reported to be in 
the a6-plane at low temperature [8]. Results of magneti¬ 
zation measurement with the field applied approximately 
along the easy axis of magnetization are shown in Fig¬ 
ure 2b. For these measurements, finely ground powders 
were mixed in epoxy and cured at room temperature in 
a magnetic field. The spin reorientation is apparent near 
170 K in Fe5SiB2. As expected, no similar anomaly is 
seen in Fe5PB2, which remains uniaxial at low temper¬ 
atures [6, 9]. In addition, the Mn and Co substitutions 
studied here appear to suppress the spin reorientation in 
Fe5SiB2. No lattice anomalies are apparent in Figure 2a 
at the spin reorientation temperature in Fe5SiB2. 


B. Magnetization measurements 

High temperature magnetization measurements were 
conducted to determine the Curie temperature (Tc) of 
each material. The results are shown in Figure 4. The 
measurements were performed in a relatively low applied 
field of 0.01 T to obtain sharp changes in J(T) at Tc- The 
Tc values are listed on Figure 4 and in Table I, and were 
determined by extrapolating to J = 0 the steepest part of 
the magnetization curve. Relatively small contributions 
to the magnetization arising from ferromagnetic impu¬ 
rity phases are seen in Fe5PB2 (Fig. 4a), Fe4CoSiB2, 
and Fe4MnSiB2 (Fig. 4b). The identities of the ad¬ 
ditional ferromagnetic phases is not clear; however, the 
magnetic impurity in Fe4CoSiB2 with Tc near 725 K is 
likely associated with the unidentified reflections in the 
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FIG. 4. Magnetization vs temperatnre measured in an applied 
field of 0.01 T for the P-containing compounds (a) and the 
Si-containing compounds (b). The data were collected on 
cooling and are normalized to the values measured at 400 K. 
Curie temperatures estimated from the data are noted on the 
figure. 


x-ray diffraction data for this material (Fig. le) based on 
diffraction and magnetization measurements on samples 
of varying purity. 

The Curie temperature of Fe5SiB2 could not be deter¬ 
mined because it is above 780 K, the maximum temper¬ 
ature for these measurements. The downturn seen at the 
highest temperatures indicates Tc is near but larger than 
this value. This is consistent with the literature report 
of Tc = 784 K for Fe5SiB2, although some change in Tc 
may be expected to occur with variation in Si/B or P/B 
stoichiometry in these compounds [6, 7]. Measurements 
on Fe5PB2 (Fig. 4a) give Tc = 640 K, similar to the 
reported values of 615—639 K [5] and 655 K [9]. 

Substitution of 20% Mn for Fe has a small effect on the 
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FIG. 5. Magnetization vs applied magnetic field measured 
at 300 K for the P-containing compounds (a) and the Si- 
containing compounds (b). Magnetization normalized by the 
saturation value is shown in (c) to facilitate comparisons of 
anisotropy fields. Note J = noM, and divide J in Tesla by 
dvr X 10“^ to convert to M in emu/cm®, or by 47r x 10“^ to 
convert to M in A/m. 


determined Tc values. In Fe 5 PB 2 Tc is increased by 10 
K, while in Fe 5 SiB 2 it is decreased by 14 K. These vari¬ 
ations are within the limits that might be expected for 
P-B or Si-B antisite defects, based on the the report for 
Fe 5 Pi_xB 2 +x [5]. A more dramatic change is observed 
for Co substitution. Replacing one of five Fe atoms with 
Co reduces Tq by more than 100 K in both Fe 5 PB 2 and 
Fe5SiB2. 

Results of magnetization measurements at 300K are 
shown in Figure 5. All materials show ferromagnetic be¬ 
havior. The saturation magnetization (Js), defined here 
as the value measured at an applied field of 5 T, for each 
material is listed in Table I. Determinations of Js were 


also performed at 5 K and the results are also listed. A 
low temperature extrapolation value of 1.73 /xs/Fe has 
been reported for the magnetic moment in Fe 5 PB 2 [5]. 
Recently Lamichhane et al. reported the saturation mag¬ 
netization of this compound to be 1.72 /xs/Fe at 2 K 
and about 910 kA/m (1.1 T) at 300 K [9]. The present 
measurements (Table I) are in reasonable agreement with 
these reports. In general, the magnetization is higher in 
the Si-based compounds than in the analogous P-based 
compounds. This is consistent with room temperature 
Mossbauer spectroscopy studies which find larger hyper- 
fine fields at the Fe nuclei in Fe 5 SiB 2 (17—23 T) than in 
Fe 5 PS 2 (15-21 T) [ 6 , 8 ]. 

Substitution of Co strongly reduces the magnetic mo¬ 
ment, as seen in Figure 5 and Table I. Comparing Js at 
5 K, a reduction of 16% and 20% are observed when 1/5 
of the Fe is replaced by Co in Fe 5 PB 2 and Fe 5 SiB 2 , re¬ 
spectively. Assuming that the average magnetic moment 
on Fe remains the same when Co is substituted gives an 
average ordered moment per Co of 0.31 fiB in Fe 4 CoPB 2 
and 0.00 in Fe 4 CoSiB 2 . The absence of a large or¬ 
dered moment on Co is consistent with our preliminary 
measurements, which suggest the absence of ferromag¬ 
netism in C 05 PB 2 . Substitution of Mn has a smaller ef¬ 
fect on the magnetization. In both Fe 5 PB 2 and Fe 5 SiB 2 , 
substituting 1/5 of the Fe with Mn reduces Js by 6 % at 5 
K. This suggests that Mn does have an ordered moment 
in these compounds. Indeed, the pure Mn analogues are 
known to be ferromagnetic [10, 11]. Assuming no change 
in the average moment per Fe atom occurs with Mn sub¬ 
stitution, the data in Table I gives 0.85 /xb and 1.27 fiB 
per Mn at 5 K in Fe 4 MnPB 2 and Fe 4 MnSiB 2 , respec¬ 
tively. The conclusion that substituted Mn is ferromag¬ 
netic while substituted Co is not may explain the sharp 
decrease in Tq that accompanies the addition of Co but 
not Mn (Fig. 4 and Table I). 

The rate at which the magnetization approaches satu¬ 
ration (Fig. 5) is determined by the magnetic anisotropy. 
This can be parameterized by the anisotropy field Ha, 
which is the applied field required to saturate the ma¬ 
terial in the magnetically hard direction. Since a poly¬ 
crystalline material with randomly oriented crystallites 
contains some fraction with the hard axis along the mea¬ 
surement direction, J ys H will also saturate near the 
anisotropy field, but the approach is typically quite grad¬ 
ual and Ha can be difficult to define precisely in this way. 
In some cases, the singular point detection technique 
[14, 15] can be used to identify anomalies in d^J/dH^ 
which occur at Ha, but that was not successful using the 
data for the present samples. This technique requires 
that there is a sharp (ideally discontinuous) change in 
the slope of J{H) at Ha in the hard direction, so that 
the second derivative contains a local minimum at this 
field. But this minimum occurs on top of a smoothly 
changing background arising from the portions of the 
sample with the hard axis aligned in other directions. 
Better estimates can typically be obtained from measure¬ 
ments on single crystals or magnetically aligned powders. 
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This was recently done for Fe5PB2 single crystals [9], and 
Ki was determined to be 0.38 MJ/m^ at room temper¬ 
ature. The anisotropy field can then be estimated by 
HqHa = 2Ki/Js to be 0.83 T. An anisotropy field closer 
to 1 T might be estimated from the data shown in Fig. 
5c, and this disagreement might be due to the neglect 
of demagnetization effects in the present analysis. With 
these limitations in mind, the data in Fig. 5c suggest 
that Co and Mn substitution reduces Ha in Fe5PB2, but 
may slightly increase Ha in Fe5SiB2. Measurements on 
single crystals would be desirable to confirm this. 


C. First principles calculations 

First principles calculations of Fe5PB2 and Fe5SiB2 
were performed to attain theoretical insight and un¬ 
derstanding of the magnetic properties of these materi¬ 
als, using the plane-wave density functional theory code 
WIEN2K [16]. We used the generalized gradient approx¬ 
imation of Perdew, Burke and Ernzerhof [17]. Eor the 
P compound we used sphere radii of 1.77, 2.15 and 2.06 
Bohr for B, Ee and P and an RK^aa; of 7.0 was used. 
Here KKmax is the product of the smallest sphere radius 
and the largest plane wave vector. Eor the Si compound 
we used sphere radii of 1.73, 2.24 and 1.82 Bohr for B, 
Fe and Si, respectively, with an HKmax of 7.0. Approxi¬ 
mately 1000 k-points in the full Brillouin zone were used 
for these calculations, with spin-orbit coupling included 
(excepting the internal coordinate optimization). 

We briefly summarize the results for the P compound 
(a more detailed discussion of these calculations can be 
found in Ref. 9). We find an average saturation mo¬ 
ment of 1.79 /rs/Fe, in good agreement with our experi¬ 
mental value; this value includes small negative moments 
—0.05 p-b) on the P and B atoms and a small orbital Fe 
moment of approximately 0.03 ^s/Fe. From the calcu¬ 
lations we find uniaxial anisotropy, as in the experiment, 
with a magnetic anisotropy constant Ki of 0.46 MJ/m^, 
which yields an anisotropy field Ha of approximately 0.9 
Tesla, in reasonable agreement with the experiment. 

The situation for the Fe5SiB2 compound is significantly 
more complex and interesting. As with the P compound 
we find a ferromagnetic state with a saturation moment 
of 1.83 /is/Fe , including negative moments of-0.06 hb 
on the Si and -0.1 on the B and a small orbital mo¬ 
ment of 0.04 fj,B per Fe. This result compares well with 
the measured value of 9.06 hb per formula unit or 1.81 
/is per Fe(Table I). We present the calculated density- 
of-states in Fig. 6 below, which depicts a substantial 
exchange splitting of 2-3 eV as well as an electronic struc¬ 
ture dominated around the Fermi level by the Fe, as ex¬ 
pected in an Fe-based ferromagnetic material. 

For the calculations of magnetic anisotropy we used 
higher precision, with as many as 10,000 k-points in 
the full Brillouin zone used in the self-consistent runs. 
Unlike in the P compound [9], however, we find pla¬ 
nar anisotropy with a first anisotropy constant Ki of 



FIG. 6. The calculated density-of-states of Fe 5 SiB 2 


-0.42 MJ/m^. Although Fe5SiB2 is a uniaxial magnet 
at room temperature, a spin reorientation occurs below 
room temperature as noted above and first observed in 
Ref. 8, so that our calculation is consistent with the 
ground state T — 0 configuration. Note that the theo¬ 
retical calculations were performed using a 300 K struc¬ 
ture (with internal coordinates optimized), so that the 
observed spin-reorientation does not appear to be the re¬ 
sult of lattice expansion but may rather originate in dy¬ 
namical effects. These two scenarios were hypothesized 
[18, 19] as competing explanations for the spin reorienta¬ 
tion [20-22] in the manganese-based ferromagnet MnBi, 
which is planar below 90 K and uniaxial above. Here we 
suggest that dynamical effects from the variation in mag¬ 
netic anisotropy associated with the motion of the ions 
away from their equilibrium positions are responsible for 
the spin-reorientation. This is an effect distinct from 
both that of lattice expansion as well as that of the gen¬ 
eral reduction of magnetic moment with increasing tem¬ 
perature, which generally reduces magnetic anisotropy. 

We have also performed calculations of the magnetic 
properties (using the virtual crystal approximation) of 
the partially Mn substituted compound Fe4MnSiB2. As 
with the other compounds we find a strongly ferromag¬ 
netic ordering, with average spin moment 1.88 /is per 
transition metal, which is slightly larger than the exper¬ 
imental value of 1.50 /is per transition metal (Table I). 
It is possible that some of the alloyed Mn anti-align with 
the Fe, which would reduce the overall moment. Unlike in 
the pure compound Fe5SiB2 we find uniaxial anisotropy 
in Fe4MnSiB2. The calculated anisotropy constant Ki 
is 0.26 MJ/m^, yielding an anisotropy field of 0.47 T, in 
reasonable agreement with the experimental results from 
Fig. 5. While a microscopic explanation for the change in 
anisotropy with Mn alloying is not readily available, we 
note that similar effects have been calculated to occur in 
other Fe-containing compounds such as Fe3Sn[23]. Thus, 
the present result should not be considered anomalous or 
unusual, and may motivate further study of Mn substi¬ 
tution as a way to manipulate the magnetic anisotropy 










in other Fe compounds. 

IV. SUMMARY AND CONCLUSIONS 

We find the tetragonal CrsBa structure-type for all 
compositions and temperatures studied. Anisotropic 
thermal expansion is observed in Fe5PB2 and Fe5SiB2, 
with a relative change in a being 4—6 times larger than 
that of c between 300 and 20 K. Fe5SiB2 in fact displays 
negative thermal expansion along the c-axis over this 
temperature range. All of the compounds have the easy 
axis of magnetization along the crystallographic c-axis at 
room temperature. The saturation moment of Fe5SiB2 
is measured to be 1.81 ^e/Fe at 5 K, in good agree¬ 
ment with density functional theory calculations which 
find 1.83 /cs/Fe. Inspection of isothermal magnetiza¬ 
tion curves at 300 K suggest that Fe5SiB2 has a smaller 
magnetic anisotropy than Fe5PB2. Calculations (T = 0) 
predict planar anisotropy for Fe5SiB2, in agreement with 
experimental findings [8], of similar magnitude but oppo¬ 
site sign from that found in Fe5PB2 [9]. The substitution 
of Co and Mn does not strongly change the strength of 
the anisotropy near room temperature, but may slightly 
decrease it in Fe5PB2 and slightly increase it in Fe5SiB2. 
An anomaly in the magnetization associated with the 
rotation of the spins into the o6-plane is observed near 
170 K in Fe5SiB2. No such anomaly is seen in Fe5PB2, 
Fe4CoSiB2, or Fe4MnSiB2 down to 30 K, showing that 
the substitution of Co or Mn for Fe suppresses the spin 


reorientation. Our calculations agree with this finding, 
showing uniaxial anisotropy in Fe4MnSiB2. Substitution 
of Co is seen to strongly suppress the ferromagnetic be¬ 
havior, reducing the magnetic moment by about 20% and 
the Curie temperature by more than 100 K. Mn substitu¬ 
tion has a much smaller effect. This is consistent with Co 
being non-magnetic (which is non-intuitive) and Mn cou¬ 
pling mostly ferromagnetically with the Fe in these struc¬ 
tures. The results presented here demonstrate chemical 
tuning of the strength and anisotropy of the ferromag¬ 
netism in Fe5PB2 and Fe5SiB2, and that first principles 
calculations can reliably predict the magnetic moments 
and easy axis of magnetization. Therefore we expect that 
continued work involving closely coupled experiment and 
theory efforts may identify compositional modifications 
that would lead to enhanced permanent magnet proper¬ 
ties in this family of compounds and others. 


ACKNOWLEDGEMENTS 

Experimental work (M.A.M.) was supported by the U. 
S. Department of Energy, Office of Energy Efficiency and 
Renewable Energy, Vehicle Technologies Office, as part 
of the Propulsion Materials Program. Theoretical work 
(D.S.P.) was supported by the Critical Materials Insti¬ 
tute, an Energy Innovation Hub funded by the US De¬ 
partment of Energy, Energy Efficiency and Renewable 
Energy, Advanced Manufacturing Office. 


[1] J. M. D. Coey, IEEE Transactions on Magnetics 47, 4671 

( 2011 ). 

[2] B. Aronsson and G. Lundgren, Acta Chemica Scandinav- 
ica 13, 433 (1959). 

[3] B. Aronsson and I. Engstrom, Acta Chemica Scandinav- 
ica 14, 1403 (1960). 

[4] S. Rundqvist, Acta Chemica Scandinavica 16, 1 (1962). 

[5] A. M. Blanc, E. Fruchart, and R. Fruchart, Annales De 
Chimie France 2, 251 (1967). 

[6] L. Haggstrom, R. Wappling, T. Ericsson, Y. Andersson, 
and S. Rundqvist, Journal of Solid State Chemistry 13, 
84 (1975). 

[7] R. Wappling, T. Ericsson, L. Haggstrom, and Y. Ander¬ 
sson, J. Phys. Colloques 37, C6 (1976). 

[8] T. Ericsson, L. Haggstrom, and R. Wappling, Physica 
Scripta 17, 83 (1978). 

[9] T. N. Lamichhane, V. Taufour, S. Thimmaiahb, D. S. 
Parker, S. L. Budko, and P. C. Canfield, “A study 
of the physical properties of single crystalline fe 5 b 2 P,” 
arXiv: 1508.05629. 

[10] D. M. de Almeida, C. Bormio-Nunes, C. A. Nunes, A. A. 
Coelho, and G. G. Coelho, Journal of Magnetism and 
Magnetic Materials 321, 2578 (2009). 

[11] Z. G. Xie, D. Y. Geng, and Z. D. Zhang, Applied Physics 
Letters 97 (2010), 10.1063/1.3518064. 

[12] J. Rodriguez-Garvajal, FullProf.2k, version 5.30, March 
2012, ILL.,. 


[13] L. Lundgren, G. Tarmohamed, O. Beckman, B. Garlsson, 
and S. Rundqvist, Physica Scripta 17, 39 (1978). 

[14] G. Asti and S. Rinaldi, J. Appl. Phys. 45, 3600 (1974). 

[15] F. Bolzoni and R. Gabassi, Physica B 346-347, 524 
(2004). 

[16] P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and 
J. Luitz, WIEN2k, An Augmented Plane Wave -I- Local 
Orbitals Program for Calculating Crystal Properties (K. 
Schwarz, Tech. Univ. Wien, Austria) (2001). 

[17] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

[18] K. V. Shanavas, D. S. Parker, and D. J. Singh, Sci. Rep. 
4, 7222 (2014). 

[19] V. P. Antropov, V. N. Antonov, L. V. Bekenov, 
A. Kutepov, and G. Kotliar, Phys. Rev. B 90, 054404 
(2014). 

[20] B. W. Roberts, Phys. Rev. 104, 607 (1956). 

[21] A. F. Andersen, W. Hagl, P. Fischer, and E. Stoll, Acta 
Ghem. Scand. 21, 1543 (1967). 

[22] M. A. McGuire, H. Gao, B. G. Ghakoumakos, and B. C. 
Sales, Phys. Rev. B 90, 174425 (2014). 

[23] B. G. Sales, B. Saparov, M. A. McGuire, D. J. Singh, 
and D. S. Parker, Sci. Rep. 4, 7024 (2014). 



